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In this paper we study the effect of finite size on the two-flavour colour superconducting state. 
As well as restricting the quarks to a box, we project onto states of good baryon number and onto 
colour singlets, these being necessary restrictions on any observable "quark nuggets" . We find that 
whereas finite size alone has a significant effect for very small boxes, with the superconducting state 
often being destroyed, the effect of projection is to restore it again. The infinite-volume limit is a 
good approximation even for quite small systems. 
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I. INTRODUCTION 

It is widely believed that quarks, the basic building 
blocks of strongly interacting matter, lose their confine- 
ment at some large density or temperature, at or near 
the chiral phase transition. Many models, for example 
the instanton liquid model, suggest that the unconfined 
quarks exert an attractive force on each other. The BCS 
mechanism predicts that an arbitrarily weak attractive 
force leads to an instability of the Fermi-surface, giving 
rise to a colour superconductor This new phase of 
quark matter has been the subject of many papers in 
recent years; see Ref. for a detailed list of references. 

The interplay of the order parameter with the flavour 
symmetry (either SU(2) or SU(3), since only the Hght- 
est quarks are relevant) can generate many interesting 
phases. These include the two-flavour colour supercon- 
ductor (2SC) which pairs quarks of different flavour and 
colour and leaves one colour unpaired |^]; the colour 
flavour locked phase (CFL) which is possible in the three- 
flavour case Q ; and the crystalline colour superconduct- 
ing phase dj. In general, a robust colour superconduct- 
ing gap of the order of 100 MeV has been predicted by 
calculations which use different interactions and differ- 
ent superconducting states. The density at which the 
maximum pairing occurs is several times nuclear density. 

An important question is what happens in finite sys- 
tems. If we ever succeed in creating such a phase in the 
lab, it must occur in finite, colour-singlet, states. Equally, 
if we have a nugget of quark matter at the centre of a 
neutron star, we expect it to be colour neutral, since it 
originates from the compression of a colour neutral object 
and is surrounded by colour neutral matter. 

The techniques with which to discuss such phenomena 
have long been used in nuclear physics, where the su- 
perfluid nature of finite nuclei should not be described 
by the number-non-conserving BCS state, but by a pro- 
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jected one [|6|. The additional complication here is the 
requirement of colour neutrality, since the Cooper pairs 
in 2SC superconductivity are colour anti-triplets. Hence 
an additional projection onto a colour-singlet state is re- 
quired. In contrast the Cooper pairs in nuclear physics 
already carry zero angular momentum. 

In this work we will use the model of Ref. ||] and work 
at zero temperature. We limit our study to the simplest 
state, the 2SC. We model the effects of finite size by 
confining the quarks to a cubic box with antiperiodic 
boundary conditions. As a result boundary effects are 
not included. After projection the exact equations for the 
superconducting gap are very involved, so in this work 
we confine ourselves to a one-parameter variation. With 
these limitations, we find that the superconducting gap 
persists after, though not always before, projection in 
systems of finite size. 

Finite size effects in superconducting matter have also 
been considered in ref. j7|, where it was found that the 
stability of small lumps of CFL matter was enhanced. In 
that paper however the emphasis was on thermodynamic 
effects; shell structure and projection onto states of good 
quantum numbers were not considered. 

The paper is organised as follows: in Section II we 
describe the model and fix the notation; in Section III 
we apply the projection over baryon number to the 2SC 
and obtain an expression for the expectation value of the 
energy in this projected state; in Section |^ we apply 
the projection over the colour singlet to the 2SC and 
generalise the expression given in the previous section. 
In Section M we give the numerical results, and finally, 
in Section W\ we state our conclusions. The reader will 
find technical details and some useful formulae in the 
Appendices ^ and 



II. THE MODEL 

We wish to investigate the effects of finite size, definite 
baryon number and colour neutrality on the colour su- 
perconducting state. As mentioned in the introduction, 
a variety of models have been so far used to study the 
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emergence of colour superconductivity. We follow ref. 
and use a standard NJL model with massless quarks car- 
rying two flavours and three colours. We leave for future 
work the study of more complex scenarios. 

Following ref. ||] we model the attractive interaction 
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^^int = 



where flavour and colour indices are Latin and Greek 
respectively, and the tensor 'Bki^ap-yS is given by 

^kl-a0jS = Efci {SSajSpS — SasSpj) . (2) 

At zero density the quarks acquire a dynamical mass, 
through the generation of a chiral condensate. This mass 
drops as the density increases, reaching zero at the chiral 
phase transition. At higher densities the colour conden- 
sate which signals the onset of colour superconductivity 
forms. 

In order to regularise the contributions stemming from 
large momenta, the authors of Ref. ||] have supplemented 
the interaction ([l]) with a form factor, to be associated 
with each quark of a given momentum. We choose to use 
a form factor of the form 

The parameter e determines the sharpness of the cut- 
off, and has been taken to be 10. This leaves two free pa- 
rameters, A and the coupling constant K in the interac- 
tion. These are fixed at zero density by requiring the pion 
decay constant to have its physical value, and by choosing 
the dynamically generated quark mass to be 400 MeV; 
this gives A = 700 MeV and K = 1.755 x 10-°^ MeV'^. 

The modes of the fields are quantised in the standard 
way, and we define the creation operators o,]^i^ ji-^ia^k) , 

j^^^^(fc), R)iai^) f*-"" particles, antiparticles and 
holes relative to the Fermi sea, Ikp), which can thus be 
written as 

IfcF) = n ^(^^^ - fc)«L(fe)«La(fc)|0) ■ (4) 

cxik 

We choose to work with states in which the Fermi mo- 
menta kp for all three colours are equal. In fact the 
properties of the BCS state are independent of the initial 
Fermi momentum for the condensing quarks. However 
colour projection becomes involved ii hp is different for 
the spectator quarks. 

The Fermi sea is the natural choice of ground state at 
low to intermediate density. At large densities the inter- 
action (|l|) is responsible for the condensation of Cooper 



responsible for the pairing of quarks at large densities 
with a contact four-quark interaction whose form is sug- 
gested by the instanton-induced interaction. The Hamil- 
tonian then has the form H = Hq + iJjnt with 
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pairs of quarks, and we get a new form for the wave func- 
tion, 

I*) = Gl Gr \kp), (5) 

where 

Gl = 

n [cos0^(fc) -f 6„^3sin^^i(fc)e'«i«a[lJfc)a[2^(-fc)" 

{aP)k 

n [cos0i(fc) -f e„03sin0i(fc)e»««WSt,,Jfc)St,2^(-fe)" 

{al3)k 

n [cos0«(fc) + e„^3sin0g(fc)e'«S(fe)ctj^Jfc)ct^2^(-fe)" 

(a/3)fc 

= Gla Glb Glc (6) 

and similarly for Gr, with the substitutions L <r^ R. 
The sum denoted (a/J) is taken over a, /? = 1,2 and 2, 1 
only; the BCS angles O^^^^ and S^i^^ c)^k) are to 
be determined variationally. 

The BCS wave function \9) of Eq. (||) describes the 
condensation of pairs of quarks with the same helicity, 
but antisymmetric in flavour and in colour, as can be seen 
by inspection of Eq. (j^) . The condensate carries a colour- 
3 order parameter. Only two colours, by convention 1 and 
2 (or red and green), condense, leaving the third colour 
inactive. This corresponds to a spontaneous symmetry 
breaking. 

The BCS angles are obtained by minimizing the ther- 
modynamic potential 

F = /i7V|^'). (7) 

The phases f are fixed by requiring that the expectation 
value of TJint is maximally attractive; this gives 

= -^i = -^B = 

^{A,B,C)'^ ^{A,B,C) = ^- (8) 

Minimising with respect to the functions 9 one finds the 
gap A and the angles d{A,B,c)if^) fo'" particles, antipar- 
ticles and holes respectively to be given by the solutions 
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-K Ekl-aPjS j (fxil)ji-^^{x)ll)Lk-i{x)lpR2p{x)ll^Lu{x) -|-h.C., 
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FIG. 1: Variational solution for the infinite system at fi — 
500 MeV. The angles are in units of 7r/4. The solid, dotted 
and dashed lines correspond to particle, antiparticle and hole 
contributions respectively. 



of the coupled equations 

tan20(A.B,C){k) 



k — fi + ^ — k 



and 



l^^— k^dkF'^(k) 

Jo 



eik-kp) sm29A{k) + 



sm2eB{k) +9{kF - k) sm2dc{k) 



(10) 



In Fig. |l] we display the solutions for the BCS angles 
Oa, Ob and 9c in the infinite volume for a density corre- 
sponding to a chemical potential ji = kp = 500 MeV. 
The solid, dotted and dashed lines correspond to the 
particle, antiparticle and hole contributions respectively. 
The solutions quickly fall off at momenta larger than the 
cutoff A as a result of the action of the form factor. Pair- 
ing is most effective at the Fermi surface, fc « /i. 

The BCS vifave function (||) docs not have a definite 
baryon number (the condensate is a superposition of 
states with different number of pairs) and it is not in a 
colour singlet state. However, the number of particles is 
fixed on average by the Lagrange multiplier /i in Eq. (|^). 
Because the fluctuations in the baryon number behave as 
^fN , they become negligible as a fraction of the baryon 
number when the condensate contains a large number of 
particles. 

In Fig. 1^ we show the quark density as a function of the 
chemical potential for a Fermi gas of quarks (solid line) 
and for the two-flavour BCS condensate (dashed line). 
At small density the condensate contains more particle 
pairs than antiparticle or hole pairs; as the density grows, 
the reverse is true, since the particle states start to be 
suppressed from the form factor and more (Fermi) hole 
states become available. 

In this paper we are interested in the effect of finite size 
on the colour superconductor. To do this we consider the 
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FIG. 2: Quark density in the Fermi gas (solid) and in the BCS 
condensate (dashed) as a function of the chemical potential. 
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FIG. 3: Solid line: number of red and green particles in the 
Fermi function of the chemical potential for a box of 

size L = 6 fm. Dashed line: expectation value of the number 
of red and green particles in the BCS condensate as a function 
of ^. At the intersections between the two curves the number 
of pairs of particles in the BCS state equals on average the 
number of pairs of antiparticles and holes. 



quarks to be in a cubic box with antiperiodic boundary 
conditions. We start with a Fermi sea in which all levels 
with k < fi are occupied, and consider the projected and 
unprojected BCS states built on this. For simplicity we 
only consider filled shells. 

In Fig. H we have plotted the quark density as a func- 
tion of fi for the two active colours in the Fermi gas (con- 
tinuous line) and in the BCS condensate (dashed line), 
in a finite box of size L = 6 fm. Notice that the dashed 
curve does not depend on kp; in fact only quantities in- 
volving the third (inactive) colour will depend onkp. Be- 
cause we work with equal values of kp (that is the same 
filled shells) for all three colours, at the intersection be- 
tween the two curves the box will contain on average the 
same number of red, green and blue particles. These are 
the points at which we will project out colour singlets. 

In the following sections we study the occurrence of 
colour superconductivity in systems of finite size and 
with a fixed number of particles. We first project the 
BCS wave function onto states of definite baryon number 
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(Section [II) and then onto colour-singlet states (Section 



III. 



PROJECTION ONTO DEFINITE BARYON 
NUMBER 



where N is the number operator, N is the number of 
particles in the Fermi sea and C is a complex number 
on the unit circle with phase (f). The rotated state can 
be created directly from the Fermi sea by the rotated 

operators G^^MO = '^G^l,r)C-^ 



There are various ways to solve the problem of picking 
a component of fixed particle number from a BCS state. 
The most efficient ones are based on the calculation of a 
contour integral; we follow here the method of residues 
of Rcf. 1^. This will be used to project the colour super- 
conducting state considered in Eq. (||) onto a state with 
a fixed number of particles. 

We start by defining the rotated BCS state 



1^(0) = eHiv-iv)0/.|^^ 



_ ^»(JV-A')0/2| 
^(JV-JV)/2|^^ 



\^{0) = GL{OGBiO\k, 



(12) 



(11) where 



Gl(C) = n [cos^^i(fc) + Cea/33sin0i(fc)e^^5('=)aLJfc)a[2^(-fc) 



n [cos0i(fc) -Hrea/33sin0i|(fc)e^«-(^)6t,,Jfc)6t^2^(-fe) 



(a/3)fe 



n [cos0«(fc) +C6„^3sin0g(fc)e*«o(fc)4jfe)cl,2^(-fc) 



(c,fe) n ^ 



(13) 



{af})k 



{al3)k 



and similarly for Gfl(C). Because the creation of antipar- 
ticles or holes lowers the baryon number, the correspond- 
ing terms in G contain C* or, equivalently, 1/C. 

The number-projected wave function is produced from 
an appropriate superposition of rotated states: 



^l*(C)) 



(14) 



where n is the number of pairs. 

The only operators whose expectation values in the 
number-projected state are non-vanishing are those 
which are particle-number conserving and so commute 
with the number operator N. This enables us to write 



(*„|0|*„) = 



(15) 



A convenient way to calculate these matrix elements 
is to find the Thouless operator S{() which maps |^) = 
into |*(C)) §: 



^(01*) = l*(C)>. 



(16) 



The explicit expression for S is given in Appendix 
There we also give the definitions of the quasi-particle op- 
erators a(C, k), /?(C, k), 7(C, k) which annihilate the state 

1^(0). 

The Thouless operator can be expanded in a series 
of terms containing zero, two, four, . . . quasiparticle 
creation operators: 



(17) 



ra=0 



Depending on the form of the interaction, only a limited 
number of these terms will contribute to the matrix el- 
ements of the Hamiltonian. Luckily, the interaction (|l|) 
selects only the first two terms in the expansion. For 
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left-handed particles these are 



Hamiltonian in the number projected state: 



qiO) 



c(l) 



(C) 



(C) 



W (cos26ii(fc) + Csin26ii(fc)) 

{a/3)fe 

4"](c) 



'LA\ 

xE 

X 



(C~ 1) £^^3 sing^(fc)cosgjt(fc)e'^^(fe) 
cos26'^(fc) + Csin26l^(A:) 

4i^(l,fc)d[2,(l,fc) . (18) 



The function S'(o)(C) = (*|*(C)) 

can be expanded as a 
Taylor-Laurent series about C = 0; with real coefficients 
dn- These are related to the normalisation constant C„ 
in Eq. (14), by 2niCn — l/v^m and to the expansion of 
the BCS state |^') in terms of states of definite particle 
number: 



I*) 



(19) 



Since |^) has unit norm, the sum of the dn is unity. 
From Eq. (nsl) we obtain the expectation value of the 



J 



(*„|i?|*„) = -27rz \Cn\ 



^(vI'|i?|vI;(C)) . (20) 



Explicitly, we obtain for the non-interacting part of the 
Hamiltonian: 

k 

{e{kF - k) [dn -f 2cOs2f?c(fc) J,,„(0c(fc))] 
+ 2e{k-kF) SIU^Oa (fc) la.n {0 A (k) ) 

+ 2smHBik)IbA0Bik))} . (21) 

In this expression the first term describes the hole con- 
tribution, the second the particle one, and the last the 
antiparticle. The first term (d„) in the square brackets 
describes the spectator colour; the second term reduces to 
twice this in the absence of pairing (when 0(^a,b,c) = 0)- 
The interaction Hamiltonian leads to the more compli- 
cated form 



(*„|7?int|*„) 

8K 



k k' 



X <j e{kF - k)0{kF - k') sm^^sm29c{k) sin 20c (fc') JbA^cik), Ocik')) 
2'^'='^' sin 20c (fc) sin20B(fc') JbA^c{k),OB{k')) 



+ 2B(kF - fc)cos 



2e{kF -k)e{k' ^kF)cos^-Y- sm26c{k)sm2eA(k') JcA^c{k),6A{k')) 



sm 



2 

2 '''='=' sin20B(fc)sin20B(fc') Jb,„(0B(fc), 0b(A:')) 



2e{k' - kF) sin^^ sin 26*5 (fc) sin 261^ (fc') JcniOsik), ^^(fc')) 
e{k - kF)9{k' - kp) sin^"^ sin 29 A{k) sin26'^(fc') Ja,ni^Aik),9Aik'))} 



(22) 



r 



The volume of the box is = L^. The contour integrals 
are included in the definitions of 2(a.b,c),n and J(a,b,c),n 
(see Appendix A); dk^i is the angle between k and k' . 
The form factors in the expression ( p^ ) suppress the con- 
tributions coming from large momenta. 

The expectation value of the number operator is ob- 
tained from Eq. ( pl[ ) by dropping the single-particle en- 
ergy k which appears immediately after the summation 
sign and swapping the sign of the antiparticle [B) term. 



These expressions (21) and ( p2[ ) reduce to the expecta- 
tion values in the unprojected BCS state if the integrals X 



and and the constant dn are set to one. The infinite- 
volume gap equations (^) and (^ are easily obtained 
by minimising {H — fj,N) with respect to the angles 9, 
if the sums are replaced by integrals. However in the fi- 
nite volume projected case, the equations are much more 
complicated. For a start, in the cubic box there is no 
reason to expect the angles to depend only on the length 
and not on the direction of k. This introduces a depen- 
dence on the specific geometry that one might be quite 
happy to ignore. However even then the dependence of 
the T's and J^'s on the 0's renders the equations very 



6 



complicated. We do not attempt to solve these. Instead 
we assume that the form of the 9s is given as before by 
Eq. (||), and minimise the expression for {H — ^N) ob- 
tained from Eqs (|2|) and (|2^) with respect to the gap A 
only. 

This process of restricting the BCS angles to have the 
same form as in the unprojected state will presumably 
introduce fewest errors where the projected state has a 
high overlap with the unprojected state. Thus although 
in principle we can project a state of a given baryon num- 
ber from any BCS state, in practice we choose values of 
/i where the expectation value of the number of pairs in 
the unprojected BCS state has the value desired in the 
projected state (see Fig. ||). 

In the BCS state there are non-vanishing expectation 
values for diquark condensates such as ipRii{x)ipii22{x). 
However as this operator changes the quark number by 
two units, it does not commute with the number operator 
N and the expectation value in the number-projected 
state will vanish. None-the-less, the square of the diquark 
condensate operators do not vanish, and they act as order 
parameters in this case. 



IV. COLOUR PROJECTION 

We now turn to the issue of colour projection. As 
we have pointed out before, the colour-superconducting 
state which forms at large densities must be a colour 
singlet. However, the state in Eq. does not fulfill 
this requirement, given the form of the operators G(i_fl). 
For this reason we want to project the colour-singlet state 
from Eq. (|^). This requires us to integrate over the group 
manifold: 



(23) 



Here dflg is the volume element on the group manifold 
and l^fra) is the number projected BCS state. Ug is a 
unitary operator which performs a rotation in the SU (3) 
space. 

Following 1^ we parametrize an element g of SU(3) in 
the form 

(24) 

where Xi are the Gell-Mann matrices, and write the vol- 
ume element as 



1 



d'^g = 7^8 8111(06) sin(0'J sin(0) sin 



28 



X d<j)a dcjib d(/)c dip d(f>^ dcf),^ dcf)^ d<j)e ■ (25) 



To obtain the realization of these rotation operators when 
acting on our quarks states, we replace the A's with the 
transition operators Qa — fd'^x'ilj^{x)^ip{x), which are 
associated with the SU(3) group. 



Because of the residual colour symmetry, the operators 
associated with the red-green SU(2) subgroup annihilate 
the BCS state, and also commute with the number pro- 
jection operation. Thus 



Q 



1,2,3 



*„) = 0. 



(26) 



We can simplify the algebra considerably by noting 
that 



Qs = I d^xij\x)^4^{x) 



Nil +N22-2 N: 



33 



(27) 



where 1, 2 and 3 are the values taken by the the colour in- 
dex. Since we are working with states built on Fermi seas 
which have equal numbers of red, green and blue quarks, 
we see that Qg annihilates only the number-projected 
BCS state |\E'o) which contains zero net pairs {i.e. the 
number of pairs of holes and antiparticles is the same as 
the number of pairs of particles). We therefore restrict 
our projection to this state. 

Thus the only operator which acts non-trivially is Q5, 
which is non-diagonal in the colour space and connects 
one of the two paired colours with the remaining one. As 
a result we are able to reduce the integration over the 
eight different angles (jji in Eq. (|2^) down to an integra- 
tion over the single angle (jj for the rotation generated by 
Q5. The residual volume element is therefore 



dfl^ = ^ sm{(j)) sin^ ( | 



= sin —Id sin — 



(28) 



where < < tt is the angle associated with Q^. 

As before, we want to project states with zero pairs 
from BCS states which have zero pairs on average, so 
we are restricted to working only at certain values of /i. 
These are the values at which the two curves of Fig. ^ 
cross. For small boxes, this is a significant restriction on 
the densities we can consider. 

Just as in the case of number projection particle- 
number-conserving operators have non-zero expectation 
values, so when we project onto a colour singlet only 
colour-singlet operators will have non-vanishing expec- 
tation values. Such operators commute with the colour 
rotation, and so we can write (cf. Eq. (|l5|)) 



(*o|0|*o) = 



JdQg {^0\Ug\^0) 



(29) 



We can write the colour-number projected state in 
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FIG. 4: Approach of the gap (left) and energy per quark for the paired colours only (right) to the infinite volume result (dashed 
line) at a fixed chemical potential fi = 500 MeV. 



terms of the colour rotated operators G'(L,ij)(C)j 

dC 



where G(i_fl,)(C,0) = UgG(^L^R){C)U},; for example, the 
operators for the left-handed particles read 



|*o> - Co 
= Co 



C 

dC 

c 



dr!5Gi(C,0)Gfl(C»|fcF) , (30) 



GLA12{k, C, < 
CLA2l{k, C, < 



cos^^(fc) + C sin^^i(fc)e'^-('=) ( cos | d[,,{k) - sin | dl,,{k) ) d[,,{-k) 



cose'i(fc)-C sin6li(fc)e*«5('=) 412^^) (cos J a[2i(-'^) - sin 



"L23 



i-k) 



(31) 



(recall that the first numerical index on the creation oper- 
ators is the flavour index, and the second is colour). This 
clearly reduce to the original form Eq. (|l^) for = 0. By 
looking at the expressions ( pl| ) we see that, by acting with 
Clai2 on the Fermi sea, one can create pairs of particles 
in the colours 12 or in the colours 23. The probability 
of each of these events will depend on the magnitude of 
the angle 0; clearly, the two extrema of the interval (i.e. 
(p = and 4> — tt) correspond to pure pairing in the 12 
and 23 colours respectively. 

Once again, calculations are most conveniently done 
with the aid of a Thouless operator, W , which fulfills the 
relation 
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(32) 



The expectation value (E9f) therefore reads 



1 



{0) = -2ni\CoY / dn^idc- (vi'|0T4^(c,</.)|*(C)). 



C 



(33) 



The explicit form of the operator W is given in Appendix 



We can expand T4^(o)(C,0) as a Laurent series in C, 
with coefficients dn{(j)). The normalization constants C„ 
are related to integrated coefficients dn = J dil^dn{4>) 
as in the pure number-projection case the C„ are to the 
Laurent coefficients o?„ of S^-^\(^). 

The forms of the expectation values of Hp, Hjnt and N 
are as given in the previous section, in Eqs (|2l|) and (|2^) , 
but setting n — 0, and with the replacement of dn with 
dn and of the integrals I and with the new integrals 
/ and J which are given in Appendix Once again 
a one-parameter variation with respect to the gap A is 
performed. 



V. RESULTS 

In this section we present our final results for the pro- 
jected states. 

First however we look at the effect of finite size on the 
unprojected BCS state. As the size of the box grows the 
fractional fluctuations in the baryon number will drop as 
1 / VN, and we expect that the value of the gap and of the 
energy should converge toward their values in the infinite 
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FIG. 5: Left plot: distribution of the Laurent coefficients of 
S^°^ (C) for a box of size L = 6 fm and for different densities 
corresponding to the intersections between the two curves in 
the previous figure. The number of particles corresponding 
to each curve is indicated. Right plot: distribution of the 
rescaled Laurent coefficients d'„ — dn/do, as functions of the 
rescaled variable n' = {do/\/2n) n. 



volume. This behaviour is evident in Fig. |j, which shows 
the superconducting gap and total energy per particle. 
Convergence is obtained quickly and continuously. 

Now we turn to the effects of number and colour pro- 
jection. As we have explained in Section IV, the pro- 



jection over colour is carried out on BCS states having 
the same total baryon number as the underlying Fermi 
gas they are built on. This restricts us to certain values 
of n, as we saw from Fig. ^. In Eq. ( p^ we introduced 
the function 5^"^ (C) whose Laurent coefficients dn are re- 
lated to the amount of the projected state with n pairs 
contained in the unprojected state (see Eq. (|l9|)). 

Fig. H shows the distribution of the Laurent coefficients 
for the function S''*'-'(C) at these points. Because we are 
working at our chosen values of fj,, the coefficients peak 
at n = 0, and we see that the unprojected BCS state has 
a substantial overlap with the projected state with zero 
pairs. The distributions are easily approximated with 
Gaussian functions. 

In Fig. ^ we show the superconducting gap as a func- 
tion of chemical potential for a range of box sizes. The 
dashed line is the finite- volume unprojected gap, and for 
comparison the continuous line shows the result for the 
infinite-volume case. We remind the reader that because 
we perform only a limited variation after projection, as 
explained above, the projection can only be done for cer- 
tain discrete values of fj, for which the average net number 
of pairs in the unprojected BCS state is zero. At these 
points we show results for the unprojected BCS (circles, 
which must coincide with the dashed line), and for the 
number-projected (triangle up) and combined number- 
colour projected BCS (triangle down). 

The most striking finite- volume effect is the disappear- 
ance of the gap for values of ^ which are far from shells. 



These regions are confined to very low fj, for boxes above 
3 fm in length. In these regions the effects of projection 
on the gap can be striking, as can be seen even at 7 fm. 
However in all other regions, where the unprojected gap 
does not disappear, projection has very limited effect. In 
general the results for colour and number projection lie 
above those for number projection alone; the latter show 
almost no difference from the unprojected case unless the 
unprojected gap vanishes. 

The effect of finite size on the energy of the condensate 
has already been shown in Fig. ^. The effect of projec- 
tion is to lower the energy by less than 1 MeV even for 
the smallest boxes considered. In the approach to the 
infinite- volume limit, shell effects for the unpaired quarks 
vanish more slowly than finite-size or projection effects 
on the condensate. 



VI. CONCLUSIONS 

In this work we have studied the effects of finite size 
over the colour superconducting state, by considering 
quark matter at finite density and zero temperature en- 
closed in a cubic box. We have considered the 2SC su- 
perconducting state of ref. ^ and we have calculated 
the effects of the projection over definite baryon number 
and over a colour singlet state. The most significant fi- 
nite size effect is the vanishing of the gap for values of 
the chemical potential between widely spaced shells. For 
these cases, projection has a significant effect, restoring 
a sizeable gap. Otherwise the effects of projection are 
slight, though in all cases each projection increases the 
gap to some degree. 

This behaviour points to the interesting possibility 
that the interplay between the chiral and the supercon- 
ducting phase can also change in a finite system. This 
hypothesis can be verified by extending the present cal- 
culation, allowing a competition between the two phases. 
This topic is currently under investigation jiot . 

There are other aspects which require future investiga- 
tion as well: it would be interesting to extend our results 
to more general form of pairing interactions, rather than 
the simple "instanton-motivated" interaction of Eq.(|]). 
For technical reasons, the number projection which was 
carried out exactly in the present work, would require 
more effort in a general case, where the interaction contri- 
butions will contain not only terms which mix the helic- 
ity, such as -0^ ipR, but also terms which connect fermions 
with the same helicity, such as ipj^ ipL- It would be more 
appropriate in such cases to apply approximate number 
projection techniques, which have been developed in the 
contest of nuclear theory |^. 

Another aspect of interest, which we plan to investigate 
in the future, is the effect of the colour projection on 
different superconducting states, and in particular on the 
colour-flavour locked phase (CFL), which was originally 
studied in . Because in this state the colour and flavour 
degrees of freedom are "locked" , the colour projection 
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FIG. 6: Gap for the finite box at a fixed size of L = 3, 3.5, 4, 5, 7, 8 fm: unprojected (circles) and dashed line, colour projected 
(cross), number projected (triangle up) and number-colour projected (triangle down). The solid line is the gap for the infinite 
volume unprojected BCS state. 



will affect both colour and flavour, and could possibly 
yield larger effects than the one observed in the case of 
the 2SC. 

We also remark that in the present analysis we have 
not performed a full variation after projection. Instead 
we have followed a simpler approach, assuming for the 
solution the same functional form as in the unprojected 
case and obtaining the numerical solution by performing 
a minimization of the total energy with respect to the 



gap A. 
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APPENDIX A: NUMBER PROJECTION: SOME USEFUL FORMULAS 

1. Quasi- particle operators 

We give here the exphcit expressions for the quasi-particle operators annihilating the state |i/'(C)) of Eq. (jl^). After 
some algebra we find that these can be written as follows: 

dLia(C,fe) = cos6'i(fc)aLia(fe)-Cea/33sin6'^(fc)e*«i('=) a[2/3(-fe), 

dL2/3(C,fc) = cos6'i(fc)aL2/3(-fc) + Cea/33sine'i(fc)e*«S('=) a[i„(A;), 

/3ma(C,fc) = cos^i^(fc)6fli„(fc)-Cea/33sin^;«(fc)e'«SW fet^^/^l-^), 

/3i?2/3(C,fe) = cos0f(fc)6K2/3(-fc)+C*ea/33sin0f(fc)e^«-W ^UW' 

7ma(C,fe) = cos0^(fc)cfl.i„(fc)-Cea/33sin0«(fc)e*«S« fiU(-fc), 

7i?2/3(C,fc) = cosd§{k)cR20{-k) +Cec.p3smeg.{k)e'iS(k) c^^^jk). (Al) 

Notice that these operators fulfill canonical anticommutation relations. The inverse equations, which allow us to 
express the original operators in terms of the quasi-particle ones, are also useful and read: 

aLia{k) = cos6'i(fc)dLi„(C,fc) + Cea/33sin6'i(fc)e^^5('=) ci[2;3(C,fc), 
aL2i3{-k) = cos6li(fc)ciL2/3(C,fc)-Cea/33sin6'i(fc)e*«^('=) d[i„(C,fc), 

bma.{k) = cos0f(fc)/3ma(C,fc) + C*ea/33sin0f(fc)e'««('=) /3]^2/3(C,fc), 
bB2p{-k) = cos0i(fc)/3fl2MC,fe)-rea;33sin0i(A:)e'«2W ^t^^^(^,fc), 
cma(fc) = cos0g(fc)7fli„(C,fc) + C*ea/33sin0g(fc)e^«S(fe) ^t^2^(^,fc), 

Cfl2/3(-fc) = cos0g(fc)7fl,2;3(C,fe)-Cea/33sin^g(fc)e'«§W 7]jiJC,fc). (A2) 



2. Thouless operator 

We now consider the operator defined in Eq. (|l6|). For the part corresponding to G\^{C,k) we look for solutions 
satisfying 



These have the forms 



k)_ 


^LAq/3(1' k) 


— ^^Aafii^^k), 


k) 


^LSq/3(1' k) 


— '-'LBa/siC^k), 


k) 


^^CapO'^k) 





(A3) 



<{C) = 1, 



<(C) = 1, 
(C) - 1, 



^«2^(C) = (C-l)tan0^(fcy«-W, 
^«2''(C) = (C-l)tan0i^(fcy«-W, 
^^2''(C) = (r-l)tan0g(fc)e^«S«. 



(A4) 



By simple substitutions wc can also obtain the similar solutions for the part corresponding to G\^{C,,k). 
Finally the Thouless operator is: 



S{0 = Sl{C)Sr{0 = Sla{C)Srb{C) Src{C)Sra{OSlb{C)Slc{0 , (A5) 
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where 



Sla{0 = 


n 

{a(3)k 


"i + (C- 


- l)e„^3tanei(fc)e'S-W a{,Jk)a{^^{- 


fe)" 


SRBiO = 


n 


"i + (C* 


- 1) e^0sta.ne^{k)e^^-^''%,^{kpj,,^{- 


-k) 


Src{0 = 


TT 

{al3)k 


1 + (C* 


- I)ea83 tan6'g(A:)e'^"('')cl,i„(fc)ct fl( 


-k) 


Sra{0 = 


n 

{a0)k 


"i + (C- 


- l)e„^3tane:«(fc)e'«-(*=) 


-k) 


Slb{0 = 


n 

{ap)k 


"i + (C 


- 1) e„^3tan^^(fc)e*«-('=)6l,Jfe)SL^(- 


-k) 


Slc iO = 


n 

{a0)k 




- l)e„^3 tan0^(fc)e^«cWc[^Jfc)c[2^(- 


-k) 



To apply Wick's theorem it is preferable to express the equations above directly in terms of the quasi-particle 
operators: 



Sla{Q 



Sra{0 



Srb{Q 



^sin' 



(a/3>fe 

^«5WdLji,fe)42/3(i,fe)' 



i*)= n [(cosxw+c 

(a/3>fe 

+ (C - l)ea/33sin6»i(fc)cos6'i(fc)e*^'*^' 

{al3)k 



{al3)k 



, (C-l)ea;33sin02(fc)cose5(A:)e 
l*)= n [(cos^^I(fc) + C*sm 

(af3)k 



I*), 



(al3)k 

= n [(cos^^£(fc)+c 

{a/3)fc 



+ (C - I)eo/33sin0g(fc)cos0f (fc)e^«S«^t^^Ji,fe)42^(l,fe) 

sin26i£(A;)) 



l*>, 



+ (r-l)eo;33sin^^|^(fc)cos0|^(A:)e^«««/3iiJl,fc)/3[2;3(l,fe)] I*), 



' sin20g(fc)) 



■^iJclC) !*>= n [(cos'^g(fc)+C* 

(a/3)fe 

+ (C - l)e„^3sin^«(fc)cos^^(fc)e*«§W7L„(l,fc)7L/3(l'fc) 



VT> -/ -n^^j ov / ov / 

n [(cos2^^(fc) + Csin20^(fc)) 

(a/3)fe 

+ (C - l)e„;33sin0^(fc)cos^^(fc)e'«SW7L„(l,fe)7L^(l,fe)" 



(A7) 
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m 



3. Integrals 

We have seen in Section [11 that the matrix elements of operators in number projected states can be expressed 
terms of a hmited number of integrals, namely 

JaAOM ^ ^ / ^^^"^(0 ^^^^20^ + Csin2^^lKcos202 + Csin^^s) ' 
JbAdiM ^ / ^S^""^ (0 ^^^^2^^ ^ sin^^i) (cos202 + C* sin2^2) ' 

Zm J 4"^^^ (cos'^t'i + 4* sm &i)(cos^t'2 + C^m 6^2) 

The integrals for the single particle operators are obtained from those above: 

IiaAc)A^) = ^(a,6,c),„(^,0) (A9) 

APPENDIX B: COLOUR PROJECTION 
1. Thouless operator 

We can now give the explicit expression for the Thouless operator which can be derived in a similar fashion 
to S. Following the same approach as before we first look for an operator, T, that is at most bilinear in the quark 
operators and that it can be applied to G{k, () to give the colour rotated operator G{k, (f). 

For example, in order to obtain the operators in Eq. (|3^) we need 

iLAi2{k,(j)) = 1 + wi a[i;^(fc)aLii(fc) + W2 a|^i3(fc)aLii(fc) 

iLA2l(fc» = 1 + W3 a^2l(~'^)"i2l(-fc) + W4 ai23(~'^)^i2l(-fc) • (Bl) 
tLA12(fc,0)G'LA12(fc,C)|fcF) = GLA12{k.X,4')\kF) 

iLA2l{k,4>)GLA2l{k,0\kF) = GLA2l{kXA)\kF)- (B2) 



The conditions 



yield the solutions 

(j) (j) 
wi = W3 = cos ~ 1; W2 — Wi = — sin — . (B3) 

Finally one can write the total operator T in the factorised form: 

T{<t>) = Y\. iLAap{k,(l>)iLBap{k,(l))iLCap{k,(l)) iiiAap{k,(l))iuBa0{k,(j>)iB.Cal3{k,(l)) . (B4) 
{a/3)fe 

As a result, the Thouless operator for the combined colour and number projection will be obtained by the application 
of f and S: 

w^(C,0) = r(0)^(c) 

= Y]_ WLAal3{k)wLBap{k)wLCap{k)wB.Aap{k)wB.Bap{k)wRCal3{k) , (B5) 

We will simply state the results for the component corresponding to left handed particles: 

Wla{c^,0 = wj°l{ct>,0 + wi}l{4,X) + --. (B6) 

where 

wi°i(0,C) = n (cos2^i(fc) + Ccos^sin2^^(fc)) (B7) 



13 



and 



E — ^7 — 2 — 71 — \ ""^(i' fe)42^(i, fc) 

^ (cos20^(fc) + Ccos|sin20^(fc)j 



^ C 81116*^ (fc) sin |e*«^W / t . ^ t . ^ t . x t . A I 



The last term pairs the third colour with one of the other two. This term does not to contribute to the expectation 
value of the Hamiltonian. 



we now 



2. Integrals 

define the following objects: 



In analogy with Sec. A3 

2Tn J C + J (cos^^i + sur6'iCcos f)(cos2 



Ccos| 



-sin^^aCcosI)' 



Jb.7l 

J, 



2^z/C"+W (cos20i+sin20iC*cos|)(cos202+sin202C*cos|) 
2mJC^+^J ^^'^^(cos20i+sin20iCcos|)(cos2(?2+sin202Ccos|) 



(B9) 

S J ^cos-(7i -f Sin (711, COS -j ^i^cos-t72 -T »m U2^ cos ^ ; 

where df^s is given by Eq. (28). These reduce to the integrals J7 of Eq. (AJ) if is set to zero and the colour-group 
integral is suppressed. 

The integrals for single particle operators will be obtained from the ones above as a particular case: 

I{a,b,c},n{9) = J(a,b,c},n{9,0). (BIO) 

In order to calculate these integrals, we write the Laurent series for (p) as 



n 

Numerically it is found that the Laurent coefficients as a function of colour angle can 
single parameter a: 



This approximation is particularly useful in the numeri cal c alculation, sinc( 
analytically. With this simplification, the integrals in Eq. (^^ ) can be written 



(Bll) 

be fit by Gaussians with a 

d„(0) « d„(0)e-t^'"'^ . (B12) 

the numeri cal c alculation, since it allows us to calculate the integrals 
;rals in Eq. (^^ ) can be written in a closed form as: 



in 



^ (sin^fli sin202A„,m+2(^?l, ^2) - cos^^i cos2^^2A„,™(^l, 02)) 



7n 



^ (sin^fli sin202A„,m(ej, ^2) " cos^^i cos2^2A«,™+2(0J , ^2)) 



m 



^ (cos^^i sin2^^2A„,„,+2(^?^, 02) + sin^^i cos202 A„,,„(0J , ^2) 



2sin20i sin^e^ A„,,„+i(0^ , ^2))] , 



(B13) 
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where 6'^ = ^tt — 9 and we have defined 



and 



A„,„i(^?i, 6*2) 

h = 
h = 



J_IJ( I 

27ri / C""™ (cos26li + sin26iiC)(cos26l2 + sin26'2C) ' 

J_IJC I 

27ri / C"-™ (cos^^i + sin26'iC)2(cos26'2 + sin^^aO^ 



.2-^ -^sin=^* 8- (a2 + 4a + 8) 



(B14) 



(B15) 



The parameter a is defined in Eq. (B12). The integrals of Eq. (B14) can be also evaluated analytically. 
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